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INTRODUCTION 


in  this  report  we  present  the  wort  performed  under  this 
contract  not  described  in  the  two  progress  reports  already  submitted. 

The  object  of  the  wort  was  to  invesyigate  factors  afecting  the  natural 
fate  of  aerosol  and  chemical  agents  during  transport  in  Western  Europe. 
We  have  investigated  rainout  occult  and  dry  deposition. 

In  section  2  of  the  report  we  describe  measurements  of  the  turbulent 
deposition  of  cloud  droplets  isize;liimi  and  dissolved  chemical  agents  to 
a  grass  covered  hill  side.  In  section  3  these  results  are  used  to 
predict  the  deposition  of  chemical  agents  to  elevated  topography  of  a 
wide  range  of  scales.  These  are  compared  with  deposition  rates 
associated  with  rain-out  using  a  model  of  the  seeder-feeder  effect. 

In  section  4  we  describe  results  of  an  experiment  to  measure  the 
turbulent  deposition  of  sub-micron  aerosol  by  eddy  correlation.  In 
section  5  we  present  the  results  of  the  UXIST  contribution  to  a  joint 
aerosol  transport  experiment  conducted  by  UMIST  and  the  University  of 
Galway.  In  section  6  the  main  conclusions  of  the  work  are  summarised. 
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Report  on  the  Measurements  of  the  Size  Dependence  of  Cl  Quo 
Droplet  Deposition  at  a  Pill  Sue. 


Abstract 


Measurements  have  been  made  of  the  deposition  of  cioud  droplets  and 
dissolved  chemical  species  to  a  grass  moorland  vegetation  rieasuremeni?  of 
water  fluxss  t'avs  DP°r!  mads  dy  tds  Qradtsn^  •jsipf’ 

Knollsnberg  FSS  prcbss  and  by  a  wSiQhmcj  jysimstsr  ucpuSit'iun  vsiocltics 
for  cloud  liquid  water  in  the  range  O.i  to  0.24  g  m*3  at  l  0  rn  were  close  to 
those  for  momentum.The  FSS  probes  enabled  the  determination  of  deposition 
velocities  as  a  function  of  droplet  side,  it  was  found  that  the  deposition 
velocity  increased  with  droplet  radius  from  2.5  pm  to  5.5  pm  radius  . 
Between  6.5  and  12.5  urn  radius  deposition  velocities  exceeded  those  for 
momentum  out  decreased  sharply  between  1 2.5  ano  15  5  urn 

in  periods  of  thin  cloud  (<0.l  g  m*3)  evaporation  from,  the  surface 
occurred  despite  continued  deposition. 
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1  Introduction 


It  has  beer.  Known  for  soma  time  that  the  direct  deposition  of  droplets, 
so-called  'occult  deposition',  (Nagel  1956;  Kerfoot  1963).  may  be  an 
important  pathway  for  the  deposition  of  pollutants  Relatively  few 
measurements  have  been  made.  Dcllard  and  Unsworth  (1983)  used  may 


irr.pactors  to  measure  gradients  in  liquid  water  content  from  which 
deposition  velocities  were  calculated.  Similar  measurements  were 
performed  by  Dollard  et  ai  (1983)  where  some  preliminary  measurements  of 
ion  fluxes  associated  with  the  liquid  water  fluxes  were  made,  in  this  paper 
we  present  detailed  measurements  of  droplet  flux,  and  the  associated 
chemical  deposition,  by  two  different  techmaues  and  investigate  the 
venation  of  deposition  velocity  with  drop’et  size 

A  large  collaborative  experiment  is  in  progress  at  Great  Dun  Fell  in 
Cumbria,  UK,  847  m.  as',  to  investigate  aqueous  phase  cloud  chemical 
processes  and  deposition  of  species  in  aqueous  solution  to  the  hill  surface. 
Results  from  the  cloud  chemistry  experiments  have  been  described  by 
Chandler  et  a!.  (1988)  and  results  from  measurements  of  the  altitudinal 


variation  of  wet  deposition  by  rainfall  in  Chcularton  et  <?/(1933)  and  Fowler 
et  a)  (1968)  Full  details  of  the  area  and  the  measurements  made  at  each 
site  may  be  found  in  these  papers 

!n  this  paper  we  describe  measurements  of  the  deposition  of  cloud 
droplets  to  the  moorland  vegetation  on  the  hillside  This  was  done  both  by 
using  the  gradient  technique  and  measuring  the  rate  of  increase  in  weight  of 
a  representative  section  of  turf  placed  on  a  lysimeter  Simultaneous 
measurements  of  the  chemical  composition  of  the  cloud  water  were  made 
enabling  fluxes  of  the  major  species  to  the  surface  to  be  deduced.  The 
measurements  were  made  on  the  SW.  face  of  Great  Dun  Fell  about  150  m 
below  the  hill  summit  which  rises  650  m  above  the  Eden  Valley  to  the  SW 
This  area  is  uniformly  covered  by  rough  grass  and  measurements  made  at 
this  site  using  Porton  anemometers  have  shewn  a  very  good  logarithmic  wind 
profile  in  the  lowest  3  m  of  the  boundary  layer  in  conditions  of  cloud  and 
strong  winds  These  measurements  indicate  that  the  site  is  well  suited  to 
the  use  of  the  gradient  technique.  Further  details  of  the  airflow  and 
turbulence  characteristics  may  be  found  in  Gallagher  et  si  (1933) 


2  Theoru 
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C'oud  droplets  ore  transported  to  the  surface  by  eddy  diffusion  and 
gravitational  sedimentation.  The  droplet  flux,  F,  within  the  surface  layer  is 
given  by  e.q.  Chamberlain  (1967)  and  can  be  written 


F  =  K  (2) 


♦  Vs  (r)  C  (z,r)  =  v'd  (z,r)  C  (z,r) 


(1) 


where  C  (z,r)  is  the  concentration  at  height  z  for  a  given  size,  K  (z)  is  the 
yorti^gi  9jjQy  pif  fustvlty,,  v$  (r)  tbs  termlnsl  velocity  of  tpo  droplet  and  Vjj 

(z,r)  the  deposition  velocity.  Since  the  flux  near  the  surface  is  a  constant 
ond  ossumiriy  risutr3l  stability,  K  (2)  —  kzu^,  k  is  Von  Karmans  ccnstant,  then 

the  solution  to  (l)  becomes 


Vd  (z.r)  = 


_ Vs  (r) _ 

{  1  -  (z/zq)  'vs(r)  ' ku*  } 


(2) 


for  which  it  is  assumed  that  C(z=Zo),  where  zq  is  the  roughness  length. 
This  result  has  been  used  to  model  the  deposition  of  chemical  species  in 
cloud  droplets  to  hill  surfaces  by  Hill  et  al.  (19875. 

!f  droplet  capture  at  the  surface  is  not  rate  limiting  and  if  u*  »  V3 

(r),  then  the  concentration  profile  is  independent  of  droplet,  s^ze,  and  the 
usual  logarithmic  profile  is  obtained  for  c  The  deposition  velocity  in  this 
case  is  then  the  deposition  velocity  due  to  momentum  only  viz,  Vm  (z)  z 
ku»/ln(z-d).  Sy  measuring  3C(z,r)/dz  and  K9=(z)  we  obtain  V'd  (z,r)  using  (1). 


3  Experimental  Procedure. 

Figure  1  shows  the  experimental  arrangement  used.  A  4  m  scaffolding 
tower  was  erected  on  the  hill  side  ,  at  an  elevation  of  575  rn  above  the 
valley  floor.  This  site  has  been  described  in  detail  by  Fowler  et  et.  (1998) 
and  is  essentially  an  elevated  moorland  site  dominated  by  tussocks  of 
Eriophorum  #?and  Juncus  sp  and  providing  extensive  homogeneous  fetches, 
greater  than  300  m,  in  most  directions.  Two  independent  methods  were 
employed  to  measure  the  cloud  droplet  deposition  process  a  gradient 
technique  and  a  direct  weighing  iysimeter  method. 
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4  Gradient  Method. 


Two  Knoilenberg  forward  scattering  spectrometer  probes,  FSSP,  were 
mounted  on  the  tower,  at  0.5  m  and  3.5  rn.  above  the  surface.  These  probes 
were  set  to  size  cioud  droplets  in  the  range  1  to  15  pm  radius  with  a 
resolution  of  05  ym  radius.  Since  mean  droplet  sizes  at  this  site  are 


typically  5  to  7  ym  radius  this  setting  is  ideally  suited  for  measurements  of 
liquid  water  content  obtained  by  integrating  resulting  number  spectra.  Both 
probes  were  aspirated  at  26.5  ms-’  using  Rotron  ducted  fans.  Variations  in 
sample  volume  which  result  from  wind  ramming  and  gust  mg  were  corrected 
for  using  the  techniques  described  in  detail  by  Choularton  si  a!.  (1986). 

One  of  the  probes  had  previously  been  compared  with  several  different 
liquid  water  content  measuring  techniques.  These  included  a  Barnes  infra¬ 
red  transmissometer  operating  on  1 0.6  ym  wavelength,  Consterame  (1984) 
and  a  laser  scattering  device  Latham  si  a!.  (1984).  The  agreement  by  these 
independent  methods  was  good  with  mean  errors  of  +  10  %,  and  extrema  of  ± 
20  *.  Thus  the  error  in  absolute  liquid  water  content  using  these  devices  is 
at  worst  +10* 


in  adopting  the  gradient  technique,  the  second  probe  was  normalized 
to  the  calibrated-  probe.  This  was  dene  by  running  both  probes  side  by  side  at 
0.5  m  above  the  surface  in  conditions  with  liquid  water  contents  in  the  range 
02  to  05  gm  m-5  Commercial  nebulizers  were  also  used  for 
intercomparison  purposes  .  The  probes  were  sampled  at  1  Hz  and  integrated 
over  typically  15  to  30  minutes  to  reduce  the  sampling  errors  in  the  larger 
si2e  ranges  of  12.5-15.5  ym  where  number  concentrations  were  low  This 


enabled  channel  by  channel  normalization  factors  to  be  obtained  The 
standard  errors  in  these  factors  were  at  most  4  %  for  sizes  up  to  10  ym  and 


15  %  for  sizes  up  to  1 5.5  ym.  Linder  low  liquid  water  content  regimes  (  <  0.1 
g  m-3)  these  errors  increased  significantly.  The  intercomparison  between 
the  two  probes  is  discussed  in  detail  later 


The  large  separation  of  the  probes  meant  that  observed  gradients 
would  be  substantially  greater  than  the  corresponding  error  in  the  gradient 
measurements  except  in  the  case  of  the  largest  droplets.  Thus  detailed 
measurements  of  gradients  in  number  concentration  as  a  function  of  drop 
size  could  be  obtained  with  confidence.  The  disadvantage  with  this  large 
separation  ,  however,  is  that  a  larger  uninterrupted  fetch  is  required  with 
additional  increases  in  advection  errors  occurring.  A  iarge  number  of 
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measurements  at  this  site  have  confirmed,  however,  that  the  wind  profile  is 
closely  logarithmic  up  to  the  level  of  the  mgner  prone 

5  Direct  vVciGhmu  Method. 


A  direct  weighing  lysimeter  was  installed  next  to  the  scaffold  tower  as 
shown  in  Figure  !.  The  instrument  and  the  logging  system  has  Seen  described 
in  detail  by  Calder  el  al.  (1SS4).  Essentially  the  system  consisted  cf  a 
Sartorius  3806MP  electronic  balance  with  a  1.0  qm  accuracy  in  30  kg.  The 
oalance  was  placed  in  a  hole  1.5  m  deep  a  plastic  weighing  tray  was  then 
placed  on  top  of  this  flush  with  the  surface  of  the  ground.  The  turf  which  had 
been  removed  prior  to  digging  the  hole  was  replaced  in  the  tray.  The 
collecting  area  of  the  tray  was  0.64 m2  Polystyrene  sections  were  mounted 
between  the  balance  and  the  tray  to  exclude  any  neat,  transfer  from  the 


balance  mechanism  to  the  tray  producing  unwanted  evaporation  effects  The 
balance  was  interrogated  by  a  microcomputer,  with  a  sampling  rate  of  1  Hz  . 
This  data  was  then  averaged  to  1  minute  before  being  written  to  magnetic 
di  sc 


S  Additional  Measurements 


Measurements  of  0,  were  obtained  from  a  set  of  precision  Porton  cup 

anemometers  mounted  on  a  3  m  tower.  The  anemometers  had  previously  been 
calibrated  in  a  wind  tunnel.  Additional  wind  speed  measurements  were  made 
at  3  m  and  10  m  using  standard  Vector  instruments  cup  anemometers.  The 
output  from  the  profile  anemometers  was  integrated  over  15  minute  periods. 

The  data  from  the  two  FSSP's,  3  and  10  m  wind  speeds  and  directions, 
temperature  at  2  m,  and  a  variety  of  air  and  cloud  water  chemical 
parameters  e  g.  pH,  aqueous  phase  H2O2,  gas  phase  Su2  and  gas  phase  O3 
were  logged  continuously  at  1  Hz  as  described  by  Choularton  et  el  (!9S8).  in 
addition  bulk  samples  of  cloud  water  were  collected  on  an  event  basis  for 
later  analysis  by  ion-chromatography. 

?  intercompartson  of  FS5  Probes  and  Treatment  of  Data. 

To  obtain  high  quality  data  in  an  experiment  of  this  kind  precautions  have 
t.o  be  token  with  FSS  probes  Unless  the  optical  alignment  of  these  probes  is 
carefully  checked  systematic  differences  between  individual  probes  will 
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occur,  i f i  oroer  to  obtain  accurate  relative  differences  oetween  two  probes 
sampling  different  regions  of  a  cloud  which  nave  undergone  different  level? 

d9pl9tlQP  D92T  tfts  Qr0UPd  Consistent  °Cr,'>r,:5!sS*2*'0r^  *5—0 rs  PVJSt 
Oulainau.  For  measurements  of  total  si  Quid  watsr  content  where  the  bu»k  o» 
the  if  quid  water  is  contained  within  a  narrow  spectrum  of  drop  sizes  the 
criteria  for  accurate  normalisation  factors  win  not  oe  as  rigorous  as  for 
comparison  of  individual  si29  classss  Ths  accuracy  o^  tftece  ^orrectr ~  ° 
T  actors  will  necessor'ily  determine  the  accuracLj  With  which  the 
concentration  gradient  can  be  measured  as  a  function  of  droplet  size 

Care  must  be  taken  to  assess  tne  various  errors  mvo'ved  m  channel  on 
chann 9 1  c c m p a rl  s c n  Thee*  a.r9 

1)  Sample  error  of  ths  FSSP. 

2)  Wind  ramming  and  gust  effects  on  the  sample  volume  at  different 
levels,  and 

3}  Errors  In  ths  normalisation  factors 

The  sample  error  can  be  reduced  by  using  sufficiently  long  averaging 
times.  Figure  2  shows  a  typical  cloud  droplet  spectrum  observed  at  this  site 
props  2  mounted  at  35  m  Th9  data  were  averaQsd  ov^r  15  minutes  "^e 
liquid  water  content,!.,  was  0.24  g  nr*  ,  the  number  concentration  .Nr  130 
cm"3,  and  the  rriean  radius  R  r  5.9  um.  The  sample  errors  for  this  run  are  uf 
the  order  0  3  %  rising  to  2.7  %  for  the  largerst  drops  Because  of  the 
depletion  in  liQuid  water  content  the  lower  probe  produced  sample  errors  of 
the  order  of  2.1  %  rising  to  15.7  %  for  the  largest  drops  when  the 
normalisation  factors  are  taken  into  account.  For  periods  where  the  liquid 
water  content  was  small  (<  0  l  g  rrrs )  the  errors  for  the  large  drops  became 
excessively  large  (>50  %).  Such  periods  were  characterised  by  variable  cloud 
upse  with  large  scale  . nhomogeneities  in  the  sampled  cloud,  in  extreme 
cases  this  inhomogeneity  could  be  observed  on  scales  less  than  the  probe 
separation  Data  from  such  periods  has  been  excluded  from  the  analysis 

pr-Qcgrt^oft  hgrg 

The  error  due  to  wind  ramming  is  corrected  for  through  the  use  of  the 
following  relation ; 


sa  =  S0(  1  ♦  — ) 


(3) 


where  Sa  is  the  actual  sample  volume,  5$  is  the  sample  volume  under 
zero  wind  conditions,  u  is  the  local  wind  speed  measured  at  the  probe  level. 
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arid  oa  is  the  aspiration  speed  of  the  fan  which  is  26.5  rn  s'!  Provided  the 
nrooes  are  orient.ea  into  wind  and  tne  wino  speed  at  eacn  levei  is  t-nown  to  a 
reasonable  ijsQree  of  accuracy  (derivpd  from  tbe  profile  dais)  then  tbe 
sample  volume  error  slioutu  be  smialt.  The  sample  volume  correction  can  also 
be  made  using  the  dropiet  transit  time  method  described  by  Chouiarton  at  at. 

0986).  Agreement,  between  the  two  methods  for  Probe  2  was  good  agreeing 
tc  within  !5  %  . 

Figure  3  shows  the  results  obtained  from  seven  10  minute 
intercomparison  runs  prior  to  the  experiment  wnere  the  probes  were  run  side 
by  side  R  is  the  mean  of  tne  ratios  of  the  concentrations  for  Rrooe  t  to 
Probs  2  as  a  function  of  radius  Ths  vsrtical  bars  rsnr^s9nt  ths  standard 
srrors  for  sach  FSS  proi)6  bin  sizs.  Thsse  vvsrs  of  ths  erdsr  01  z.  tor  drops 
up  to  iO  urn  and  rose  to  typicaiiy  i5  %  for  15um  crops  When  the  prooes  are 
separated  vertically  the  data  from  eacn  channel  of  probe  l  is  multiplied  by 
the  appropriate  normalisation  factor.  The  liquid  'water  gradient,  or  the 
gradient  in  the  number  concentration,  in  each  channel  between  the  two 
heights  C8n  be  calculated.  The  errors  being  generated  by  combining  the 
sample  errors  from  the  two  probes  and  tme  error  in  R 


3  experimental  Results 


Table  1  shows  details  of  several  runs  obtained  during  the  early  nours  of 
the  morning  of  i  i  m  June  1937  over  a  +hree  hour  period  The  wind  during  this 
period  was  approximately  10  ms-’,  from  030°.  Temperatures,  measured  at 
2.0  m,  varied  from  3  8  to  4.0 «  C.  Figure  4  is  a  trace  of  liquid  water  content 
calculated  from  the  droplet  histograms  recorded  by  the  FSSP  mounted  at  3.5 
m  for  the  period  00:00  to  12:00  hrs  3.M.T.  The  large  variation  in  cloud  base  is 
evident.  Mean  liquid  water  contents  varied  from  near  zero  to  240  mqm'S 
Values  of  u»,Zo  and  d  were  obtained  from  the  Porton  anemometer  tower.  In 


each  case  very  good  logarithmic  profiles  were  obtained  as  shewn  e.g.  in 
Figure  5.  Unfortunately  during  the  experiment  two  of  the  profile 
anemometers  ceased  to  operate  and  information  was  only  available  from 
four  levels.  A  large  number  of  similar  measurements  at  this  site  with  the 
full  complement  of  anemometers  show  excellent  logarithmic  profiles 
consistent  with  the  measurements  presented  here.  Table  2  contains  the  bulk 
deposition  parameters  calculated  for  each  run  it  was  found  that,  for  periods 


of  low  liquid  water  content  e.g  Runs  2  end  8  it  became  difficult,  to  determine 
a  reliable  profile  of  concentration  due  to  very  large  short  term  fluctuations 


7 


Ki  liquid  water  content.  This  is  rsT'ecteu  bu  tne  ratios  of  v?  to  wrn  which  are 
much  greater  than  )  \f  these  two  runs  are  igr.nreo  then  t/>p  ratio 
v*/Vrn  r  0.96  i  0.12. 


figure  S  shows  the  measursmarits  obtained  i  r jTn  LiiS  il 


JIIUCiCI 


penod  00:34  to  02:00  .  Frequent  calibration  of  the  instrument  was  made, 
particularly  when  the  site  went  out  of  cloud  During  trns  mittgi  period  the 
increasing  trend  due  to  water  transfer  from  the  cloud  to  the  balance  is 
clear,  however,  the  noise  on  the  trace  due  to  the  variability  in  the  cloud 
meant  that  the  error  on  the  mean  slope  was  large  for  individual  runs,  as  is 
indicated  in  Table  3.  where  because  of  evaporation  a  negative  trend  could  be 
oDssrved  ip.  som 9  cssss  Th.9  srrors  shewn  for  th9  ^ysimstsr  srs  *h9  s * ond2rr 
errors  of  ths  slops. 

Better  agreement  could  be  obtained  for  periods  of  cloud  where  conditions 
were  relatively  steady  e  g  combining  runs  3  &  4  gave  a  lysimeter  flux  of  1 1 
±  3  mg  m~2  £-1  compared  to  a  cloud  droplet  flux  of  8.8  t  1.1  mg  nr  2  $-! 
i his  corrasponded  to  a  period  of  substantial  liguid  water  content  which 
remained  reasonably  constant  so  that  little  evaporation  from  the  surface  of 
the  lysimeter  occurred.  During  this  penod  no  precipitation  was  reported 
Even  very  light  dri zzl9  will  domip.sts  the  wster  flux  to  th9  lysimeter  hut 
will  not  Qiisct  the  gradients  from  the  FSci  probes  e  g.  xun  c; .  Figure  6  disc 
shows  1  minute  averages  of  iiauid  water  content  obtained  from  probe  2  for 
the  first  five  runs  The  slow  response  is  evident,  for  instances  when  large 
changes  in  liquid  water  content  occur  These  effects  were  particularly 
evident  for  runs  7  to  10  as  a  result  of  lew  li  quid  water  content  and  pcssiuic 
evaporation  effects 

9  Variation  of  Deposition  velocitu  with  Droplet  Size 


Each  si2e  class  was  treated  in  a  similar  manner  to  the  histogrammei 
liquid  water  content  date  from  each  probe,  using  the  channel  by  channel 
normalisation  values  shown  in  Figure  3.  Table  4  shows  the  results  from  rus 
3.  This  run  was  typical  for  periods  where  the  liquid  water  content  at  1  hi 
was  between  0.13  and  0  17  g  nr?  ie  runs  3,4,5  &  iO  To  compare  thr 
venation  from  run  to  run  the  deposition  velocity  was  normalised  using  th, 
momentum  deposition  velocity  as  shown  in  Figure  7. 

The  relationship  between  deposition  velocity  and  drop  size  showe, 
consistent  behaviour  from  run  to  run.  The  deposition  velocity  (Vg)  increased 
with  particle  sice  within  the  size  range  2  5  to  6  5  pm,  from  about  tO  mm  s 
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1C 


to  50  mm  s“ •  .  At  about  6.5  urn  droplet  radius  the  deposition  veiociiu  is 
approximately  equal  to  the  aeoosition  velocity  for  momentum  Vrn  (Figure  7). 


The  small  droplets  (2.5-5. 0  y,m)  are  not  captured  at  the  surfaces  of  the 
vegetation  as  efficiently  as  momentum,  which  is  transf erred  by  skin  friction 
and  form  drag  in  droplet  deposition  velocity,  ana  the  gradual  increase 
reflects  the  increasing  efficiency  of  impaction  processes  as  the  mass  of  the 
droplet  increases,  CftamCerlain  (1975).  The  observed  deposition  velocity 
increases  to  values  greater  than  Vm  for  the  size  range  5.5  to  10.5  urn  radius 
after  which  a  sharp  decline  is  observed  tor  1 0.5  to  15.5  urn  drooiets.  For  the 


S’ze  range  12.6  pm  to  1 5.5  pm.  droplet  deposition  velocities  are  sma'ler 
tft3n  Vrfi,  3] though  the  error  h^rc,  fnr  t hp$p  rjrcp!ets  ?re  'srje 


10  Chemical  Deposition 


Samples  of  the  cloud  water  were  collected  and  stored  in  a  ^e^,,~’n9rQtor 
at  4  0  C  for  later  analysis  Dy  ion-chromatography,  using  the  1  i quid  water 
fluxes  from  Table  2  the  chemical  deposition  rate  for  the  C'emical  agent  was 
calculated  This  is  shown  i r  Table  5 


i  \  Discussion 


The  experiments  have  shown  that  at  Great  Dun  Fell  the  deposition  velocity 
of  cloud  liquid  water  to  the  hillside  is  ulose  to  the  value  of  momentum.  This 
is  in  broad  agreement  with  the  results  of  Dullard  ft  si.  (1983)  The  measured 
chemical  agent  deposition  rate  was  i  mg  m'*  nr  .  This  may  be  compared 
with  values  typically  3-6  mg  m’*  hr  *’  reported  by  Chculartcn  gt  si.  (1969) 
due  to  rainfall  when  the  seeder-feeder  mechanism  is  operating.  The  direct 
deposition  of  cloud  droplets  is  unlikely  to  contribute  a  major  fraction  of  the 
water  deposited  to  a  hillside  The  mean  annual  wind  speed,  computed  over  a 


10  year  period,  for  Great  Dun  Fell  is  1 1  m  s' 


10  m  above  the  summit. 


With  226  days  of  cloud  per  year  each  of  4  hours  duration  this  would  give 
between  50  8nd  130  mm  of  water  deposited  to  the  hill  by  cloud  compared  to 
typically  2000  mm  by  rain.  However,  because  the  total  duration  of  in  cloud 
conditions  for  hills  such  as  Great  Dun  Fell  is  much  greater  than  the  duration 
of  rainfall,  the  deposition  of  chemical  species  by  turbulent  deposition  is  at 
least  comparable  to  that  for  rainfall  at  high  altitude  Between  the  period 
1972  and  1979  Great  Dun  Fell  was  in  c'oud  for  parts  of  226  days  per  year 
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compared  to  153  rain  daus.  This  would  result  in  cloud  oeoosit’.cn  duration 
typically  5  to  iO  times  tnat  for  ratotgii 


Tho  rnwi'ortrstirtf'C  of  pH  .^hQiv>»r-V  cnor'oc  ;ra  >-«ar-ar*Mi' 

l.iw  ■FVnvt.i-.i  Uv.vl.W  w.  L...  W..W....WN-.  $  w  •  •  - «  -  •  •  3 


•  r  l*i  h;nhur 

. Sr*~’ 


i «•  Ciuuu  vvfltsr  tnan  m  rain  water  Fowler  et  si  \  19335-  Trie  data  presented 
above  also  show  that  when  the  cloud  is  thin  the  iysirneter  may  show  a  net 
weight  loss  despite  a  significant  downward  flux  of  cloud  water  deduced  dy 
the  gradient  technique.  This  suggests  that  vapour  deficit  at  the  wet  surface 
is  causing  evaporation  of  the  water  on  the  affected  surfaces  and  provides 
exoerlmental  suoport  to  the  sugoestlun  made  by  Unsworth  (1953)  that 
simultaneous  turoulent  deposition  and  evaporation  o’'  water  from  the 
2f?9Ct9d  surfaces  Will  r9SUlt  arthaproii  "‘f  :"j"sp}T|i£2] 

species  on  the  surfaces  . 

Trie  data  oresented  in  this  paper  have  for  the  first  time  demonstrated 
the  size  dependence  of  tne  deposition  ve'ocitu  of  cloud  droplets  over  a 
natural  surface  T^e  ’.ncrsase  ’n  '•'e^odt11  b^tsyoop  rirr>rstot  farm  *> f 

2.5  and  6.5  m.71  may  be  attributed  to  the  reduction  of  the  surface  resistance 
in  this  range  These  results  are  consistent  with  the  model  predictions  of 
31mn  (1982)  for  solid  particles  it.  is  d,ff’cult  to  model  such  effects 
accurately  due  to  factors  such  as  sheltering  coefficients  ,  ey.  Underwood 
v»9S/‘).  Droplets  between  6.5  and  10.5  jj.m  idti'nuS  exceed  the  momentum 
deposition  velocity  due  to  increased  efficiency  of  collection  by  imoaction  to 
the  surface.  Above  10.5  pim  a  sharp  decrease  in  deposition  velocity  occurs 
for  radii  up  tc  ^5  5  p.00  (*he  s^e  Tv^n.  o**ppt  ic 

to  quantify  due  to  the  large  sample  errors  in  the  data  for  this  size  range  Put 
is  qualitatively  consistent  with  bounce  off  effects  predicted  for  solid 
particles  in  a  similar  range  range  and  obser/ed  by  Cnamberlain  (1967).  The 
data  may  suggest  that  droplets  in  this  size  range  are  not  being  efficiently 
captured  by  the  surface.  This  effect,  however,  is  not  likely  to  become 
important  as  in  this  size  range  sedimentation  velocities  are  beginning  to 
become  important  and  will  clearly  dominate  for  still  larger  droplets  More 
work  is  required  to  quantify  these  effects. 

The  results  obtained  here  suggest  that  the  model  of  Hill  et  a l  (1SS7) 
may  be  used  to  give  reasonable  estimates  of  turbulent  deposition  rates  to 
hills  They  assumed  that  the  droplet  deposition  velocity  was  equal  to  tnat 
for  momentum.  This  assumption  is  only  in  error  for  very  small  droplets  of 
•  adius  r.  5  pun.  this  size  will  be  achieved  within  100  in  of  cloud  base  on  hil>s 
in  maritime  locations  where  occult  deposition  is  common  More  detailed 
calculations  including  surface  resistance  effects  will  be  reouired  in  regions 
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frequently  sucjected  to  different  CCii  distributions  nmere  large  numcers  of 
very  small  droplets  are  present 
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Table  3  .  Deposition  velocities  calculated  from  Run  5  as  a  function  of  r>zs 
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mm  s'1 
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Figure  1 


Figure  2. 


Figure  3. 


Figure  4 


Figure  5 
Figure  6 


Legends  to  Figures 


The  experimental  arrangement  showing  .  Scaffold  tower 
withFSSP's,  Profile  anemometer  tower,  Lysimeter; 
and  the  Instrumented  Vehicle  &  Generator 

Histogram  of  number  concentration  (cm” versus 
radius  (pm)  for  a  typical  in  cloud  period  observed 
by  the  FSSP  mounted  at  3  5  m  . 

The  normalisation  factor  R  derived  from  running  both 
probes  side  by  side  with  the  mean  error  for  all  runs 
as  a  function  of  droplet  si2e 

Trace  showing  1  minute  averages  of  liquid  water 
content  calculated  from  droplet  spectra  observed  by 
FSSP  2  at  3  5  m  for  the  morning  of  1 1/6/87 

Wind  profile  measured  at  4  levels  for  Run  5  on  1 1/5/37 

The  variation  in  liquid  water  content  obtained  from  the  FSSP 
(dark  line)  and  the  mass  of  the  lysimeter  (dashed  line) 
from  00  00  to  02  00  GMT  (runs  2  to  5). 


Figure  7  Deposition  velocity  normalised  by  the  momentum 


Vt 

deposition  velocity  —  ,as  a  function  of  droplet  radius 
» m 

for  runs  3,4,5  &  1 0  on  1 1  /6/87. 
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HuDELLIHG  THE  DEPOSITION  OF  CHEMICAL  AGENTS  BY  VASHGUT  AND  TURBULENCE 
Abstract 


It  has  recently  been  shown  iron  model  predictions  and  measurements 
at  Great  Dun  Fell,  U.K,  that  the  variation  oi  chemical  agent  deposition 
with  altitude  associated  with  rainiall  is  much  larger  than  had 
previously  been  thought.  In  this  paper  the  model  oi  deposition  by  the 
seeder-feeder  mechanism  is  extended  to  cover  a  wider  range  oi 
atmospheric  conditions  and  hill  sizes  and  the  results  are  compared  with 
a  model  oi  turbulent  deposition  .  It  is  shown  that  turbulent  deposition 
rates  are  about  a  factor  of  5  to  10  less  than  the  seeder-feeder 
mechanism.  The  patterns  of  both  types  of  deposition  are  strongly 
affected  by  atmospheric  stability  and  windspeed. The  model  of  deposition 
by  washout  of  a  hill  cap  cloud  by  rain  predicts  that  for  a  steep  hill, 
characteristic  length  B  km,  the  maximum  rainfall  and  chemical  agent 
deposition  rates  are  downstream  of  the  hill  simmit.In  this  case,  the 
raindrop  trajectory  through  the  highest  liquid  water  content  region  is 
short,  but  the  condensation  rate  is  large  due  to  strong  updraughts. Hence 
the  greatest  effect  is  the  wind-drift  of  the  raindrops  and  consequently 
the  rainfall  and  deposition  maxima  are  displaced  downstream.  .  On  a  longer 
hill,  about  10  im  in  length,  cloud  water  is  not  so  rapidly  resupplied 
due  to  reduced  updraught. The  rainfall  rate  in  the  region  of  the  summit 
is  greater  than  for  a  steep  hill  because  the  region  of  high  water 
content  is  large  in  horizontal  extent  and  the  effects  of  wind-drift  are 
small. As  a  result,  the  maximum  rainfall  and  deposition  rates  occur 
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somewhat  upstream  oi  the  hill  summit,  depending  on  ataospnenc 
stability. 


1.  Introduction 


In  this  section  of  the  report  we  describe  a  model  of  the 
deposition  of  chemical  agents  to  elevated  topography  by  the  direct 
turbulent  deposition  oi  cloud  droplets  from  the  cap  cloud  to  the  surface 
and  washout  by  the  seeder-feeder  process. 

a.  Cap  Cloud  Kodel  and  turbulent  deposition 

This  model  is  a  development  of  that  described  in  Hill  et  al  1966 
and  so  it  will  only  be  outlined  here.  In  order  to  calculate  the 
chemical  agent  concentration  in  the  droplets  at  a  particular  point  on 
the  hill  side  we  require  a  representative  transit  time  from  cloud  base 
to  a  measuring  site  on  the  hill  surface. We  therefore  firstly  require  an 
airflow  model  of  the  hill. The  wind  velocities  need  to  be  predicted  close 
to  the  hill  surface  with  some  accuracy  as  well  as  throughout  the  cloud 
depth  because  all  measurements  are  made  within  10  metres  of  the 
ground. We  use  the  model  of  Carruthers  and  Choularton  19S2  which  was 
developed  specifically  for  Great  Dun  Feil  and  which  lias  received 
experimental  verification  at  this  site. 

We  solve  Helmholtz's  equation  for  lnviscid  flow 

Vi’S  = 0 


where  ju  is  the  Scorer  parameter  S. 


in  layer  i. 


and  Ui  is  the  geostrophic  wind,  s  iS  the  streamline  displacement. 


and  use  f \x)  as  the  lower  boundary  condition  where 

fix)  =  H 

+  **' L*  ) 

H  is  height  of  hill  and  L  half  width  at  half  height. 

The  buoyancy  frequency  H,  is  defined  as 

I  g  d6,\  » 

^0,d2^ 

* 

where  0i  is  the  mean  potential  temperature  of  the  layer  and  g  the 
acceleration  due  to  gravity. 

The  atmosphere  may  be  divided  up  into  several  layers  and  this  is 
important  because  the  stability  of  the  air  above  the  turbulent  boundary 
layer  can  have  a  large  effect  on  the  flow  of  air  in  this  lowest  layer 
<the  layer  in  which  the  cloud  frequently  resides). 

Close  to  the  hill  surface  the  Reynolds  stress  and  viscous  terms  cannot 
be  left  out  of  the  momentum  equation  and  therefore  in  the  region  just 
above  the  ground  the  wind  velocities  are  calculated  by  the  method 
described  in  Jackson  and  Hunt  1975. It  is  not  proposed  to  discuss  this 
model  in  detail  here.  It  has  recently  been  suggested  by  Nason  1586  that 
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|  the  Inner  layer  depth  within  which  the  turbulence  Is  in  local 

equilibrium  may  be  significantly  smaller  than  suggested  by  Jackson  and 
|  Hunt.  Measurements  of  the  wind  profile  using  a  3m  tower  of  precision 

anemometers  and  other  anemometers  at  3m  and  10m  above  the  terrain  on  the 
summit,  together  with  measurements  of  the  turbulence  using  a  Sonic 
Anemometer  have  been  made  at  Great  Dun  Fell  (Gallagher  et  al  1967)  . 
These  measurements  show  that  the  inner  layer  depth  is  about  30m  with  the 
turbulence  close  to  local  equilibrium  in  the  lowest  10m  at  least.  All  of 
the  hills  discussed  below  have  inner  layer  depths  at  least  as  large  as 
GDF.Hearly  all  the  depletion  in  cloud  droplet  concentration  occurs  well 
within  this  depth  and  so  this  has  no  significant  effect  on  the 
discussion  presented  below. 

Given  the  wind  field  over  the  hill  surface,  the  next  step  is  to 
incorporate  the  cap  cloud. In  this  model  the  cloud  either  begins  when  the 
rising  air  has  cooled  to  saturation  or  we  can  specify  cloud  base  as 
occurring  at  some  height  based  on  observation. The  depth  of  the  cloud  is 
estimated  from  soundings  in  the  vicinity  of  the  hill. The  density  changes 
due  to  release  of  latent  heat  during  the  condensation  process  are 
insufficient  to  alter  the  dynamics  for  a  hill  of  this  size. 

The  position  of  any  streamline  and  the  wind  velocity  along  any  point 
of  it  can  be  calculated  and  used  to  find  the  liquid  water  content  and 
droplet  size  distribution  variation  over  the  hill  using  the  standard 
droplet  growth  equations. Hence  we  can  produce  droplet  size  distributions 
and  droplet  composition  concentrations  for  a  given  position  on  the 
hill. We  are,  in  effect,  employing  a  lagrangian  parcel  model  at  this 
stage  <  containing  the  microphysics)  and  following  its  development  as  it 
is  driven  over  the  hill  by  the  wind-field  determined  as  described  above. 
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The  concentration  ox  chemical  species  may  be  allowed  to  vary  within 
the  collection  of  droplets  according  to  the  equations  set  out  in  section 
2.  and  gases  exchanged  between  the  droplets  and  the  surrounding  air. For 
the  purposes  of  this  paper  the  details  of  the  chemistry  associated  with 
a  particular  droplet  class  are  not  important  and  the  mass  of  eg. 
chemical  agent  deposited  is  determined  from  the  total  mass  present  at  a 
given  time. The  bulk  concentrations  of  the  various  ions  are,  however, 
determined  by  integration  of  ten  droplet  size  categories. The  droplet 
growth  equations  used  may  be  found  in  Hill  and  Choularton  1G86. 

The  boundary  layer  in  which  the  cap  cloud  resides  is  turbulent  so 
that  droplets  tend  to  mix  vertically  during  transit. The  approximate 
distance  through  which  a  droplet  will  move  in  time  at  is  given  by  U*ht 
where  U*  is  the  friction  velocity . Given  the  approximate  transit  time  of 
cloudy  air  from  cloud  base  to  the  point  of  measurement,  also  at,  we  can 
arrive  at  a  vertical  distance  such  that  droplets  throughout  this  depth 
will  have  an  equal  chance  of  being  intercepted  at  the  measuring  site. 

Figure  1  shows  schematically  one  such  region  of  cloud  that  this 
argument  defines,  stretching  back  upstream  from  a  particulir  point  on 
the  hill  towards  cloud  base. An  average  cloud  parcel  trajectory  can  now 
be  calculated  and  hence  a  vertical  wind  profile  for  use  in  determining 
droplet  growth. 

Knowing  the  composition  of  the  cloud  which  is  in  contact  with  the 
hill  allows  us  to  estimate  the  turbulent  deposition  rate  of  a  chemical 
species  to  the  hill  surface  of  known  roughness  length  at  a  given 
position,  A.Ve  solve 


K^zldCiz)  +  V«.C(.z)  =  i/di  z)C<zj  =  constant 

dz 

Where  C(z>  is  the  droplet  concentration  and  Vd'z;  their  deposition 
velocity . K (z)  is  the  eddy  diffusivity  =  KzU*,  and  the  droplet  terminal 
velocity  V«.<<  U.«  and  U*  is  defined  by 

U (z) =U*ln (z/Zo^ 

K 

K  is  the  van  (Carman  constant. 

Solving  far  Va(z)  at  a  point  A,  assuming  the  droplet  flux  is  independent 
of  height  near  the  ground,  gives 

Va(z)  =  U*« 

inCz/Zo) 

In  the  absence  of  detailed  experimental  data  to  the  contrary  we  assume 
that  the  deposition  is  limited  by  turbulent  diffusion. This  is  likely  to 
be  a  reasonable  assumption  since  the  average  size  of  the  cloud  droplets 
is  5  to  10  pm  unless  we  are  very  close  tD  the  region  of  cloud  base  near 
the  hill  surfaceThe  liquid  water  content  at  A  is  calculated  assuming  no 
depletion  (ie.  z  must  be  >  10m  if  z.=  =  0.02  m)  and  hence  we  use  the 
concentrations  predicted  by  the  model  10m  above  the  ground. If  Zo  =  0.5m 
,  corresponding  to  a  forested  hill,  then  z  =  25m,  U*  is  about  1ms  '  for 
windspeeds  of  15ms  '  employed  by  the  model  and  hence  an  insignificant 
fraction  of  time  is  spent  by  droplets  in  this  region  compared  with  the 
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transit  time  over  the  hill,  consequently  the  existence  oi  the  depletion 
zone  will  not  affect  the  chemistry  oi  the  deposited  droplets. 


3.  Results  of  the  Turbulent  Deposition  Model 

Table  t  presents  a  set  of  standard  conditions  used  as  input  to 
the  model.  The  input  parameters  are  typical  oi  values  measured  at  Great 
Dun  Fell  and  the  stability  parameters  chosen  give  a  supercritical 
airflow  regime  characteristic  of  hills  whose  height  is  comparable  to  the 
boundary  layer  depth.  The  hill  height  of  665m  chosen  corresponds  to  the 
height  of  GDF  above  local  terrain.  The  hill  half  length  L  of  2  km 
corresponds  to  GDF  and  10km  to  a  longer  hill. 

The  Great  Dun  Fell  hill  is  grass  covered  and  experiment  has  suggested 
that  a  roughness  length  of  2  cm  is  appropriate. However,  this  can  be 
changed  in  the  model  in  order  to  describe  other  terrain  eg.  a  forested 
hill. The  rate  of  deposition  will  be  higher  over  such  a  surface  if  Zo  is 
taken  as  50  cm  for  such  a  hill.  This  assumes  that  all  the  terrain  in  the 
model  domain  is  forested  so  that  turbulent  diffusion  delivers  the 
droplets  to  the  tree  tops  where  they  are  intercepted.  If  a  transition  in 
the  surface  properties  occurred,  eg  from  grass  to  forest  then  locally 
enhanced  deposition  would  occur  due  to  interception  at  the  forest  edge 
and  the  turbulence  would  be  disturbed  for  a  short  distance  downwind 
while  i’t  adjusted  to  the  higher  roughness.  These  complications  are  not 
treated  in  this  paper. 
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The  Short  Hiii 


In  figures  2 a  and  2b  curve  a  shows  the  occult  precipitation 
rate  and  chemical  agent  deposition  rates  respectively  for  a  short  h  i  i.  I  - 
It  can  be  seen  that  the  precipitation  rate  increases  toward  the  hill  top 
as  the  liquid  water  content  and  windspeed  increase.  The  airflow  revime 
is  supercritical  and  consequently  the  deposition  velocity  continues  to 
increase  on  the  lee-side  of  the  hill. This  does  not,  however,  produce  anv 
marked  shift  in  the  maximum  'occult'  precipitation  rate  away  from  the 
summit  because  the  water  content  decreases  as  the  cloud  moves  down  the 
lee-side.  The  chemical  agent  deposition  rate,  though,  continues  to 
increase  across  the  hill  as  the  windspeed  increases,  this  process  being 
independent  of  the  cloud  liquid  water  content  above  cloud  base  ‘where 
the  deposition  is  assumed  to  be  limited  by  turbulent  diffusion). 

Doubling  the  windspeed  ‘curves  o)  roughly  doubles  the  deposition  rates 
as  the  total  depletion  of  the  cloud  is  very  small.  Raising  cloud  base  to 
500ra  considerably  reduces  the  precipitation  rates  by  reducing  the  cloud 
liquid  water  content  but  does  not  affect  the  chemical  agent  deposition 
rates  except  by  reducing  the  part  of  the  hill  enveloped  by  cloud  <  the 
deposition  rates  of  sub-cloud  aerosol  are  assumed  to  be  small  due  to 
surface  resistance).  Increasing  the  chemical  agent  aerosol  loading 
proportionately  increases  the  deposition  rate,  curve  d.  Curve  e  shows 
the  effect  of  a  totally  neutrally  stable  atmosphere.  The  major  effect  of 
this  is  to  reduce  the  windspeeds  on  the  lee-side  of  the  hill  and  hence 
reduce  the  deposition  rates  there. 
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Curve  g  shows  the  effect  of  simulating  a  forested  hill  with  a  roughness 
length  of  0.5m.  This  results  in  an  increase  in  the  deposition  rates  of 
2.5  to  3  times  on  the  short  hill.  The  depletion  of  chemical  agent  and 
liquid  water  are  still  small 

In  order  to  illustrate  the  effect  of  changes  in  atmospheric 
stability  we  show  the  results  of  weakly  stable  upper  air  and  a  neutral 
atmosphere  with  symmetrical  flow  (figure  3>.It  can  be  seen  that  the 
stability  of  the  upper  layers  is  important  in  determining  the  deposition 
rate  especially  on  the  lee-side. Dif ferent  atmospheric  cona.  can 

result,  in  differences  of  a  factor  of  2. 

The  Long  Kill 

Figures  4a  and  4b  show  the  occult  precipi tat  ion  rates  and 
chemical  agent  depositions  over  the  hill  with  L=10kra. The  curves 
correspond  to  those  of  figures  2a  and  2b.  The  pattern  of  results  is  very 
similar  to  the  short  hill  case  with  some  exceptions.  The  deposition  and 
precipitation  rates  are  not  as  high  in  the  long  hill  case  for  potential 
flow  (curve  e),  particularly  near  the  hill  top  and  on  the  lee-side 
because  the  speed-up  of  the  wind  is  significantly  reduced  due  to  the 
larger  aspect  ratio  (see  Carruthers  and  Choularton  1982) . However  in 
supercritical  flow  regimes  the  pattern  is  more  complex  with  larger 
deposition  rates  well  downstream  of  the  summit. 

A  further  effect  is  that  for  the  depths  of  cap  cloud  chosen  (400  and 
800m),  and  when  the  hill  Is  forested,  with  a  roughness  length  of  0.5m,  a 
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significant  depletion  of  both  liquid  water  and  chemical  agent  occurs. 
This  means  a  reduction  in  deposition  rate  especially  downwind  of  up  to 
about  o054  and  30'A  respectively  when  this  is  allowed  for.  In  order  to 
simulate  this  depletion  we  remove  a  fraction  ,  F  , of  droplets  per  time 
interval,  at,  given  by 


F  =  V.3  ,  at 


He 


Where  He  is  the  depth  of  the  cloud. We  assume  implicitly  that  turbulent 
mixing  occurs  through  the  depth  of  the  cloud  in  a  time  short  compared 
with  the  transit  time.  In  all  other  cases  discussed  above,  the  depletion 
of  the  cap  cloud  is  small  and  so  its  vertical  depth  does  not  play  a 
significant  role. 

The  effect  of  Increasing  the  aerosol  chemical  agent  loading  is  to 
increase  the  deposition  rate  proportionately  everywhere  over  the  hill. 

4 . Washou  t  Mode  1 

The  other  process  by  which  cloud  water  is  deposited  onto  the  hill 
surface  is  rainfall . Because  the  average  droplet  lifetime  in  deep  cap 
cloud  is  typically  ten  to  twenty  minutes  there  is  insufficient  time  to 
generate  precipitation  sized  droplets  by  coalescence. Instead  any  rainout 
is  likely  to  be  due  to  the  sweeping  out  of  droplets  by  rain  falling  from 
higher  level  cloud  often  associated  with  a  front. This  process  is 
frequently  referred  to  as  the  ’seeder-feeder’  process,  Bergeron  1965. The 
rainfall  rate  will  therefore  be  enhanced  over  the  surface  of  the  hill 
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and  the  enhancement  will  be  related  to  the  path  of  the  raindrops  falling 
through  the  cloud  -  being  proportional  to  the  liquid  water  content 
integrated  along  the  path. The  model  is  based  on  that  of  Carruthers  and 
Choularton  1983. The  concentration  of  aerosol  chemical  agent  upwind  of 
the  hill  is  specified  and  assumed  to  be  independent  of  height  through 
the  depth  of  the  feeder  cloud.  This  may  sometimes  be  unrealistic 
especially  when  the  cap  cloud  extends  into  the  inversion  layer  as 
discussed  below.  Vhen  the  cloud  forms  it  is  assumed  that  all  of  the 
chemical  agent  is  incorporated  into  the  cap  cloud  droplets  by  nucleation 
scavenging. Typical  cloud  base  updraughts  for  this  model  are  about  1ms  ' 
and  under  these  circumstances  most  aerosol  size  categories  are  activated 
suggesting  that  in  this  case  this  assumption  is  valid. The  details  of  the 
washout  model  are  presented  in  Carruthers  and  Choularton(1983> .  It  has 
been  modified  here  to  include  calculations  of  the  washout  of  the 
chemical  agent  and  is  outlined  below. 

The  air  within  the  depth  of  the  feeder  cloud  is  covered  by  a 
streamline  following  grid.  If  the  air  following  a  particular  streamline 
is  initially  below  the  height  of  cloud  base  no  condensation  or  chemical 
agent  scavenging  takes  place  until  it  has  been  displaced  vertically  to 
the  height  of  cloud  base.  Condensation  along  streamlines  initially  above 
cloud  base  starts  as  soon  as  they  suffer  a  vertical  displacement.  The 
rainfall  rate  from  the  seeder  cloud  together  with  its  chemical  agent 
concentration  is  specified  and  the  raindrop  size  distribution  is  assumed 
to  follow  a  Marshall-  Palmer  distribution  (the  model  has  also  been  run 
using  a  Best  distribution  for  the  seeder  cloud  raindrop  spectrum  but 
this  does  not  affect  the  results  significantly  see  Carruthers  and 
Choularton  1983).  This  rain  falls  into  the  top  of  the  feeder  cloud. 
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Starting  at  the  most  upstream  element  of  the  top  row  of  the  grid  the 
liquid  water  content  increase  is  calculated  as  the  air  parcel  travels 
adiabatically  across  the  grid.  This  liquid  water  is  depleted  by  the 
scavenging  of  the  cloud  droplets  and  the  chemical  agent  they  contain  by 
the  raindrops, with  a  collection  efficiency  close  to  unity  as  the  air 
parcel  travels  across  the  grid  element.  Hence  the  nett  change  in  the 
liquid  water  content  and  the  depletion  of  chemical  agent  across  the  grid 
element  may  be  calculated.  This  process  is  repeated  along  the  top  row  of 
the  grid. 

The  increase  in  radius  and  chemical  agent  concentration  in 
each  raindrop  size  category  is  calculated  from  the  amount  of  cloud  water 
scavenged  as  the  raindrop  falls  vertically  through  the  grid  element.  The 
windspeed  and  terminal  velocity  of  the  droplets  are  used  to  calculate 
the  drift  of  each  drop  si2e  category  during  this  process  and  the 
droplets  are  moved  downstream  and  allocated  to  the  appropriate  grid 
element  of  the  next  row  down. 

The  whole  process  described  above  is  repeated  until  the 
raindrops  reach  the  hill  surface.  It  is  assumed  that  the  scavenging  of 
chemical  agent  is  negligible  for  liquid  water  contents  of  less  than 
0.02g  m  3.  The  results  are  insensitive  to  this  figure.  All  nucleation 
scavenging  will  take  place  rapidly  whilst  the  liquid  water  content  is  of 
this  order  and  below  this  value  droplets  will  exist  as  unactivated 
aerosol  for  which  the  collection  efficiency  by  raindrops  is  very  small. 

A  steady-state  situation  is  envisaged  in  which  cloud  is  being 
continually  swept  out  by  the  rain  and  replaced  by  fresh  cloud  due  to 
condensation  in  the  air  flowing  over  the  hill. The  path  of  the  droplet  is 
determined  by  its  fall  speed  and  the  flow  speed  of  the  air  over  the 
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hill. There  are  therefore  two  main  factors  which  determine  the  deposition 
pattern  .Firstly,  there  is  the  path  length  of  the  rain  drop,  determined 
by  the  hill  length  and  the  wind  velocity  and  secondly  the  rate  at  which 
water  depleted  by  the  washout  process  is  replaced  by  condensation. The 
condensation  rate  is  determined  by  the  updraught  speed  and  consequently 
hills  of  greater  aspect  ratio  and  hence  higher  vertical  winds  close  to 
the  summit  will  have  associated  cap  clouds  which  are  not  so  seriously 
depleted  of  water  close  to  the  summit.  Ve  can  summarize  these  arguments 
in  terms  of  two  length  scales; 

The  horizontal  drift  of  a  raindrop  radius  r  falling  at  speed  V,.  through 
a  cloud  of  depth  Zc 


U,  c  ZcUo 
W 

The  scale  length  of  horizontal  variation  in  liquid  water  content  q 

ha  ^  b 


Where  L  is  the  half  length  of  the  hill.  For  a  short  hill,  La*U>. 
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5. Results  of  the  Hadel  of  the  Seeaer-reeaer  effect 


a.  The  Rainfall  Distribution 

To  illustrate  these  effects  figures  5a  and  5b  show  the  rainfall 
distribution  patterns  for  the  case  of  a  short  hill  (applicable  to  Great 
Dun  Fell)  characteristic  length  2km  and  for  a  long  hill  ,  length  10km. 
Other  parameters  are  varied  as  for  the  occult  precipitation  discussed 
above. 

When  the  hill  Is  steep,  figure  5a,  the  trajectory  of  the  raindrops 
through  the  highest  water  content  region  is  relatively  short  but  the 
condensation  rate  is  large  due  to  the  large  vertical  winds. Hence  the 
greatest  effect  is  wind-drift  and  the  maximum  rainfall  rate  is 
displaced  downstream. For  a  supercritical  flow  regime  (curve  a>,  the 
maximum  enhancement  is  about  10%  larger,  and  slightly  nearer  the  summit, 
than  for  potential  flow  (curve  b).The  effect  of  wind-drift  is  markedly 
reduced  by  a  lower  windspeed  (curve  c),  or  a  shallower  cap  cloud  (curve 
dl,  which  also  has  the  effect  of  reducing  the  rainfall  enhancement 

Higher  windspeeds  (curve  e)  tend  to  reduce  the  maximum  enhancement, 
while  Increasing  the  distance  downstream  over  which  enhancement  occurs. A 
higher  rainfall  rate  (curve  f)  does  not  increase  the  enhancement  pro¬ 
rata  as  the  raindrop  sizes  also  Increase. 

In  the  case  of  the  longer  hill.  Figure  5b,  cloud  liquid  water 
scavenged  is  not  so  rapidly  resupplied  due  to  the  reduced  updraught . In 
the  region  of  the  summit,  the  rainfall  rate  Is  greater  than  for  a  steep 
hill  because  the  naximum  liquid  water  content  region  is  relatively 
large  in  horizontal  extent  and  the  effects  of  wind-drift  small  .  The 
maximum  enhancement  occurs  slightly  upstream  of  the  summit  for  the 
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supercri t ical  flow  regime  (curve  a>  and  near  the  summit  for  potential 
flow  (curve  bi.The  enhancement,  under  supercritical  flew  is  also  about 
15%  less  than  for  potential  flow. 

A  shallower  cap  cloud  reduces  the  enhancement  markedly  (curve  C', 
whereas  a  higher  wlndspeed  (curve  d)  increases  the  enhancement  due  to 
higher  vertical  winds  and  larger  condensation  rate  in  the  feeder  cloud. 

The  effects  on  enhancement  of  a  high  rainfall  rate  (curve  e'  is 
less  than  for  the  steep  hill  due  to  the  lower  cap  cloud  liquid  wa*-er 
content . 

b.  The  Deposition  of  Chemical  Agents 

[f  we  now  consider  the  deposition  of  the  aerosol  chemical 
agent  dissolved  in  the  cLoud  water  then  the  patterns  for  the  short  and 
long  hills  are  represented  in  figures  6  and  7.  It  is  immediately  apparent 
that  the  deposition  rates  are  5X  to  10X  higher  than  the  corresponding 
occult  deposition  rates 

For  the  short  illl,  figure  6  the  deposition  curves  peak 
downwind  of  the  summit.  The  changes  in  deposition  patterns  due  to 
varying  physical  parameters,  closely  follow  the  distribution  of  rainfall 
enhancement.  This  is  because,  in  each  case  only  a  small  proportion  of 
the  chemical  agent  in  the  cloud  is  washed  out.  Increasing  the  seeder 
cloud  precipitation  rate  to  5mm  h~'  (curve  e)  increases  the  peak 
deposition  rate  from  4 . A  mg  nr^h*'  to  17.6  mg  nr-h"'  and  slightly 
changes  the  distribution  of  the  deposition  due  to  the  larger  raindrops 
associated  with  the  higher  rainfall  rate. 


hi 


/ 

v* 
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For  the  long  hill,  figure  7,  there  is  a  faster  decline  in 
deposition  rate  downwind  of  the  hill  (eg.  curves  a > . This  is  because, 
although  the  water  content  is  maintained  to  some  extent  by  condensation 
on  the  upwind  side, the  aerosol  ehemicai  agent  is  not  replenished  and 
hence  declines  continuously  with  position  on  the  hill.  This  becomes  a 
significant  effect  on  the  long  hill  because  of  the  longer  transit  times 
of  air  parcels  through  the  feeder  cloud.  Again  It  can  be  seen  that,  in 
the  absence  of  chemical  agent  production,  increasing  the  aerosol 
chemical  agent  loading  entering  cloud  base  or  increasing  the  seeder 
rainfall  rate  (figure  7a,  curve  e)  has  the  most  dramatic  effect  on  the 
deposition  pattern.  Increasing  the  seeder  cloud  rainfall  rate  to  5mm  h-' 
increases  the  peak  deposition  rates  from  4.0  to  14.2  mg  nr-  Ir  '  , 
smaller  than  for  the  short  hill  due  to  the  higher  total  chemical  agent 
scavenged. 

A  potential  flow  regime,  (figure  7,  curve  b),  tends  to  produce 
a  more  symmetrical  deposition  pattern  while  a  higher  windspeed,  (figure 
7b,  curve  c),  increases  the  peak  deposition  slightly,  and  also  shifts 
the  peak  more  towards  the  summit. 

6. Discussion 

The  object  of  this  paper  was  to  compare  occult  and  rainfall 
deposition  rates  under  similar  conditions. Results  show  that 
precipitation  by  the  seeder-feeder  process  is  5-10  times  as  efficient  if 
typical  concentrations  of  chemical  agent  and  rainfall  rates  as  measured 
In  the  0.K  are  used. Many  hills,  however,  are  enveloped  by  cloud  for  much 
longer  periods  than  the  total  duration  of  rainfall  and  so  'occult' 
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deposition  tan  make  a  substantial  contribution  to  the  total  chemical 
agent  deposition. 

The  results  presented  above  have  shown  that  occult  precipitation 
rates  are  predicted  to  be  a  maximum  close  to  the  hill  top  but  with  a 
supercri tical  airflow  regime,  characteristic  of  hills  with  a  comparable 
height  to  the  boundary  layer  depth,  the  chemical  agent  deposition  rates 
are  a  maximum  on  the  lee  side  of  the  hill  for  both  long  and  short 
hills.  <It  should  be  noted  that  the  area  over  which  material  is  deposited 
is  much  less  than  when  the  washout  process  is  operating. ) 

A  forested  hill  has  a  deposition  rate  2.5  to  3  times  higher 
than  on  a  grass  covered  hill.  This  is  reduced  somewhat  over  a  long  hill 
where  depletion  of  the  cloud  by  the  deposition  may  be  important. 

The  pattern  of  deposition  during  rainfall  is  much  more 
strongly  dependant  on  hill  size.  On  the  short  hill  the  peak  deposition 
rate  and  rainfall  occur  downwind  of  the  hill  summit  due  to  the  important 
effects  of  wind-drift.  On  the  long  hill  the  peak  rainfall  and  deposition 
rates  occur  at  or  somewhat  upwind  of  the  hill  summit.  Maximum  deposition 
rates  tend  to  be  somewhat  smaller  than  on  the  short  bill  due  to 
depletion  of  the  feeder  cloud  by  the  scavenging,  although  the  total 
deposition  on  the  hill  is  considerably  larger. 
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Table  1 


Parameter 


Value ( lnit lal  il 

variable) 


stability  at  layer  1 
Stability  of  layer  2 
Stability  of  layer  3 
Height  of  cloud  base 
Height  of  hill.  H 
Half  width  at  half  height 
Roughness  length,  z.„ 

Geostrophic  wind,  U-a 
Temperature,  T  (at  cloud  base) 
Rainfall  rate (washout  model) 

Sulphate  concentration  in  seeder  ram 
Height  of  inversion 
Vertical  Depth  of  Feeder  Cloud 
Thicsness  of  inversion 


0 

2x10  •'  nr 
1x10  •  m  ’ 

300  m 
685  m 

2000  is  and  10000m 
0.  02  a 
15  ms~ ’ 

5-C 

1  mmh  ' 

20pK 
400  m 
800m 

400  m 


Chemical  Agent  loading  of  aerosol  enter 


ing  cap  cloud  1.5x10 
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Legends  to  Figures 


Figure  1 

Schematic  diagram  showing  turbulent  deposition  to  a  point  •  a  1  on  the 
hill.  The  atmosphere  has  been  divided  into  three  layers,  the  lowest  of 
which  contains  the  cap  cloud. —  —  represents  the  outer  boundary  of  the 

region  whose  droplets  may  be  deposited  at  *  a  1  . - -  represents  the 

average  path  to  '  a  '.Axes  in  km. 

Figure  2a 

Cloud  water  turbulent  deposition  rates,  D,  for  short  hill  case. 

Curve  a  Standard  conditions  (see  table  2) 

Curve  b  Geastrophic  windspeed  doubled  to  30ms  ' 

Curve  c  Cloud  raised  to  500a 

Curve  d  Chemical  agent  loading  increased  to  5  pgm  -isame  as  curve  a> 
Curve  e  Seutral  atmosphere 
Curve  g  Forested  hill. zo=0.  5m, 

(curve  f  not  shown,  same  as  curves  a  and  d) 

Rates  are  in  aah  1 . Xaxis  in  km. Dashed  line  indicates  extent  of  cloud 
when  cloud  base  is  500m. 


Figure  2b 

Sulphate  turbulent  deposition  rates  ,  S, for  short  hill  case. 
Curve  a  Standard  conditions  tsee  table  2) 

Curve  b  Geostrophic  windspeed  doubled  to  30ms  1 
Curve  c  Claud  raised  to  500m 

Curve  d  Chemical  Agent  loading  increased  to  5  pg®"-’ 

Curve  e  Neutral  atmosphere 
Curve  g  Forested  hill .  z.,=0. 5m. 

Rates  are  in  mgm  ~.X  axis  in  km  from  centre  of  hill. Dashed  line 
indicates  extent  of  cloud  when  cloud  base  is  500m. 


Figure  3 

Deposition  rate  of  chemical  agent,  S,  to  the  hill  surface  for  three 
different  atmospheric  conditions. 

curve  '  a  *  is  for  a  typical  supercritical  flow  using  the  input  shown  in 
table  2. 

curve  '  b  '  is  an  example  of  a  sub-critical  flow  pattern  <p-=l .  5x10"-’, 
^3=1.5x10'-’  m"  '  height  of  layer  2  Is  1200m  and  thickness  300m). 
curve  '  c  '  has  0  stability  in  all  3  layers. X  axis  in  km  from  centre  of 
hill.Y  axis  is  deposition  rate  in  mgsr-^h'1 
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Figure  4a 

Cloud  water  turbulent  deposition  rates,  D,  for  long  hill  case. 

Curve  a  Standard  conditions  L=10kffl(see  table  2) 

Curve  b  Geostrophic  windspeed  doubled  to  30ms  ‘ 

Curve  c  Cloud  raised  to  500m 
Curve  e  Neutral  atmosphere 
Curve  g  Forested  hil  1 .  z<,=0.  5m. 

Sates  are  in  mmh  ' . X  axis  in  km. Dashed  line  indicates  extent  of  cloud 
when  cloud  base  is  500m. 

Figure  4b 

Chemical  Agent  turbulent  deposition  rates,  s,  for  long  hill  case. 
Curve  a  Standard  conditions  (see  table  2> 

Curve  b  Geostrophic  windspeed  doubled  to  30ms  ’ 

Curve  c  Cloud  raised  to  500m 

Curve  d  omitted  for  clarity  (  curve  offscale) 

Curve  e  Neutral  atmosphere 
Curve  g  Forested  hi  1 1 .  z,,=0. 5m. 

Curve  g'  Forested  hill.  z,,=0.5m.  Ho  droplet  depletion. 

Rates  are  in  mgrn-.X  axis  in  km  from  centre  of  hill. Dashed  line 
indicates  extent  of  cloud  when  cloud  base  is  500m. 
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Figure  5a. 

Precipitation  rates.  P,  for  a  2kra  hill. 

Curve  a  :  standard  (dashed  line  see  table  2) 
b  ;  potential  flow 
c  :  5  ms-1  wind 
d  :  400m  thick  cap  cloud 
e  :  30  ms-1  wind 

f  :  5  mm  h-1  precipitation  rate  (x  a  scaling  factor  of  0.1  , 
dotted  line) 

Rates  are  in  nrnih -  '  m  - .  X  axis  in  km  from  hill  centre. 

Figure  5b 

Precipitation  rates,  P,  for  10km  hill 
Curve  a  :  standard  (dashed  line  see  table  2) 
b  :  potential  flow 
c  :  400m  thick  cap  cloud 
d  :  30  ms-1  wind 

e  :  5  mm  h-1  precipitation  rate<x  a  scaling  factor  of  0. 1  ,  dott 

line) 

Rates  are  in  mmh"'nr-.X  ax.s  in  km  from  hill  centre. 
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Figure  6 

Chemical  Agent  deposition  rates,  S,  by  seeder-feeder  mechanism  for  short 
hill . Physical  parameters  are  varied. 

Curve  a  :  standard  conditions  <see  table  2) 

Curve  b  :  Potential  flow,  neutral  atmosphere. 

Curve  c  :  30  ms  1  geostrophic  wind 
Curve  d  :  4 00  m  cloud  thickness 

Curve  e  :  5  mmh"'  seeder  precipitation  rate  <  x  a  scaling  factor  of  0.5) 
Rates  are  in  mgh’a“-.X  axis  in  km  from  hill  centre. 

Figure  7 

Chemical  agent  deposition,  S,  rates  by  seeder-feeder  mechanism  for  long 
hi! 1 . Physical  parameters  are  varied. 

Curve  a  :  standard  conditions  (see  table  2> 

Curve  b  :  Potential  flow,  neutral  atmosphere. 

Curve  c  :  30  ms-1  geostrophic  wind 
Curve  d  :  400  m  cloud  thickness 

curve  e  :  5  muh"’  seeder  precipitation  rate  >  x  0.5) 

Rates  are  in  mgh'nr-.X  axis  in  km  from  hill  centre. 
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Introduction 

This  experiment  formed  part  of  the  wider  scale  lnvesti gat ions  into 
the  dry  deposition  of  particles  and  atmospheric  pollutants  carried  out 
by  the  UKIST  Atmospheric  Physics  Research  Group. 

The  data  acquisition  was  carried  out  at  the  Research  Group's  field 
station  at  Great  Dun  Fell  <GDF)  in  Cumbria,  on  May  26th  1987.  The 
primary  aim  was  to  study  the  rates  of  deposition  of  aerosol  particles 
using  the  eddy  correlation  technique. 

Dry  Deposition  and  Eddy  Correlation 

Dry  deposition  is  the  removal  of  particles  or  gases  from  the 
atmosphere  through  the  delivery  of  mass  to  the  surface  by  non- 
precipitation  atmospheric  processes,  and  the  subsequent  chemical 
reaction  with,  or  physical  attachment  to,  vegetation,  soil,  or  the  built 
environment  ‘Dolske,  Gatz  1985). 

Eddy  correlation  is  the  measurement  of  the  net  turbulent  flux  at  the 
height  of  the  sensing  equipment.  It  is  assumed  that  the  sampling  takes 
place  within  the  surface  boundary  layer,  as  the  theory  requires  that  the 
flux  divergence  between  the  surface  and  the  sensors  is  small. 
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The  particle  flux  is  given  by: 

particle  flux  =  w* f  =  wj  -  w  I 

where  w  =  vertical  wind  speed,  and  I  =  particle  concentration.  The 
□verbar  denotes  the  mean  value  over  a  period  of  time,  and  '  indicates 
the  variable  component  of  the  quantity. 

However,  it  is  usual  to  assume  that  w  =  0.  This  is  not  usually  the 
case  in  practice,  often  due  to  misalignment  of  the  vertical  anemometer 
sensor.  It  Is  therefore  necessary  to  mathematically  correct  for  this 
when  processing  the  data.  Ve  can  realign  all  our  data  so  that  w  and  v 
are  both  zero  <v  is  the  transverse  wind  speed). 

Beumann  and  Hartog  (19S5)use  the  following  correction  for  the 
vertical  particle  flux: 

w' S'  =  (wB  -  w  B>  cos8  -  CuB  -  u  S)  cos8  sina 
-  (vl  -  v  S)  cose  si  no 


where  8  =  arctani  w _ l 

I  (u*  +  v1)1*! 

and  a  -  arctani  xl 

l  u  I 

and  u  =  longitudinal  wind  speed 

The  deposition  velocity  is  now  given  by: 
v.„  =  alSl 

i 

The  values  of  v„  for  aerosol  particles  are  generally  very  small  and 
it  is  therefore  all  the  more  important  to  determine  the  level  of  error 
we  have  in  our  calculations.  The  following  treatment  of  errors  is 
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suggested  by  Katen  and  Hubbe  (1985).  The  variance  of  the  covariance  of 
w'  and  S'  Is  given  by: 


=  C  (W  >3(J'  )Z1  -  (w'B'  )2 

The  first  term  on  the  right  hand  side  Is  the  expectation  value  of  <w‘)2 
and  (S')2.  The  second  tera  is  the  square  of  the  covariance.  The  standard 
error  of  an  Individual  measurement  of  w' S'  Is  given  by: 

£  <w'  S'  )  =  e.  ..... 


n* 

where  n  Is  the  number  of  observations  used  In  calculating  the 
covariance.  The  standard  error  of  the  deposition  velocity  Is  then  given 
by: 


£  (vd)  -  v...  ... 

S  n“ 

S  Is  the  average  particle  concentration  during  the  sampling  period  and 
Is  assumed  to  be  nearly  constant. 


Our  data  also  allowed  us  to  calculate  values  of  the  Honln-Obukhov 
length  L,  the  heat  flux,  the  friction  velocity  u,,  the  stability 
parameter  z/L  and  the  roughness  length  z0. 

The  Monln-Qbukhov  length  Is  given  by  (Bason  and  King,  1984): 


kgw'T' 

where  T  Is  the  absolute  temperature,  k  is  the  von  Karmann  constant 
(=  0.4),  and  g  Is  the  gravitational  acceleration  9.Skgms~2). 
w'T'  Is  the  covarlaDce  of  the  vertical  wind  speed  and  the  temperature 
and  is  given  by: 

w'T’  =  wT  -  w  T  =  heat  flux 

pCc, 
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p  is  the  air  density,  Cp  is  the  specific  heat  of  air  at  constant 
pressure.  This  quantity  also  has  to  be  rotationally  corrected  into  the 
mean  wind  direction. 

The  friction  velocity  ue  Is  given  by: 
u«a  =  -u'w' 

The  stability  parameter  is  given  by: 

z. 

L 

where  z  is  the  height  at  which  sampling  takes  place. 

The  roughness  length  zQ  is  given  by: 

Zo  *  _ Z _ 

I  ku  <z )  +  al 
expl  u,  l 

for  z  <  0  :  a  =  21nl  1  ♦  xl  +  lnl  1 .+  x2l  -  2arctan  x  +  s. 

I  I  2  I  1  2  I  2 

x  =  I  l...r..,15zi  '* 

I  Li 


for  z  >  0  :  a  =  -4. 7z 
L  L 

Th&-fixperlaefltal  £et~Pp 

An  Active  Scattering  Aerosol  Spectrometer  Probe  (ASASP)  model  300A, 
manufactured  by  Particle  Xeasurement  Systems  (PHS),  was  used  to  measure 
particle  concentrations.  The  ASASP  allows  particles  in  four  different 
size  ranges  to  be  observed.  The  range  chosen  for  this  experiment  was 
0.4pm  -  1.0pm  diameter.  This  range  is  sub-divided  into  15  categories, 
each  of  0.04pm  width.  The  ASASP  was  sampled'at  4Kz. 

The  wind  velocity  was  measured  with  a  Kaijo-Denki  USAT  310  sonic 
anemometer.  This  anemometer  provides  wind  speeds  in  three  mutually 


perpendicular  directions  i.e.  u.v  and  w.  The  anemometer  also  provides 
temperature  data.  The  sonic  anemometer  was  sampled  at  40Hz. 

The  anemometer  and  ASASP  were  mounted  on  separate  tripods.  However, 
eddy  correlation  theory  requires  that  the  wind  velocity  sensor  be  as 
close  to  the  particle  sensor  as  possible.  This  was  attempted  by 
Interlocking  the  legs  of  the  tripod.  The  centre  of  the  anemometer  was 
2.44m  above  the  ground,  0.45m  to  the  side  of,  0.40m  in  front  of,  and 
0.36m  above  the  Inlet  of  the  ASASP  (total  distance  from  Inlet  Is  0.70m). 
(A  similar  set-up  Is  shown  in  Plate  1.) 

The  data  from  the  sonic  anemometer  and  the  ASASP  were  logged  by  a 
HASCOH  data  acquisition  system,  and  then  transferred  to  magnetic  tape. 
The  data  were  processed  using  an  IBJI-compatlble  Tandon  computer,  using  a 
specifically  developed  PASCAL  program  to  compute  the  rotational 
corrections  and  covariances,  and  finally  the  values  of  particle  fluxes 
and  deposition  velocities.  (The  anemometer  data  had  first  to  be  averaged 
up  to  4  Hz.)  These  last  two  quantities  were  calculated  both  for  each 
Individual  size  category,  and  also  as  bulk  quantities  for  the  whole  size 
range. 

It  would  have  been  desirable  to  have  an  extensive  uniform  upwind 
fetch.  Formally  at  GDF  data  acquisition  takes  place  In  south-westerly 
airflows,  the  direction  la  which  the  terrain  Is  considered  to  be 
reasonably  uniform.  (This  direction  is  also,  conveniently,  perpendicular 
to  the  axis  of  the  ridge  of  which  GDF  forms  a  part.  This  Is  of 
particular  Importance  when  modelling  the  airflow  over  the  ridge.)  During 
the  experiment  the  wind  direction  was  easterly.  Icwever  it  ts  known  that 
for  the  purposes  of  eddy  correlation  this  wind  direction  Is  acceptable. 


The  Measurements  were  Initially  made  in  two  runs:  1642Z-1753Z  and 
1755Z-2102Z.  On  processing,  the  data  were  split  into  four  runs,  each  of 
an  hour  duration:  1642Z-1742Z,  1755Z-1855Z,  1656Z-1956Z  and  1957Z-2057Z. 
These  times  correspond  to  the  transition  between  late  afternoon  and 
early  evening  (including  sunset). 

The  effect  of  this  is  clearly  shown  by  the  presence  of  clear  trends 
in  the  data  (see  table  of  results).  The  heat  flux  changes  from  around 
lOOVa'-2  upwards  during  the  first  run  to  downwards  during  the 

final  run.  The  stability  changes  from  the  unstable  value  of  -0.046  in 
the  first  run  to  neutral/stable  values  in  the  other  three  runs. 

Zo  changes  significantly  during  the  period  of  the  experiment. 
Although  the  stability  changes  may  have  had  a  small  effect,  the  main 
cause  of  the  change  Is  thought  to  have  been  a  change  in  wind  direction. 

The  particle  fluxes  for  individual  sizes  vary,  often  dramatically, 
from  category  to  category.  It  was  decided  that  it  was  therefore 
desirable  to  calculate  the  bulk  fluxes  and  deposition  velocities  for  the 
whole  size  range.  The  bulk  flux  of  particles  changes  from  around  40 
partlclescc~'s~’  downwards  in  the  first  run  to  100  partlclescc* ’ s'' 
upwards  In  the  final  run.  This  trend  Is  mirrored  by  the  bulk  deposition 
velocity  which  changes  from  0.07cms"‘  downwards  in  the  first  run  to 
0.17cas~'  upwards  in  the  final  run. 

This  accumulation  of  data  trends  tends  to  suggest  that  the 
transition  from  an  afternoon  of  11°C  to  a  cool  evening  of  ?°C  has  had  a 
marked  effect  on  the  dry  deposition  of  the  aerosol  particles  under 
observation.  Ve  should  therefore  be  very  wary  of  drawing  any  solid 


conclusions  fron  our  results.  The  clearest  Indications  that  we  should 
not  place  our  total  trust  in  the  results  are  the  error  calculations 
which  suggest  an  error  In  the  results  of  up  to  an  order  of  magnitude. 

Some  of  this  error  Is  most  likely  to  have  bees  caused  by  the 
unsatisfactory  positioning  of  the  sonic  anemometer  and  the  ASASP  in 
relation  to  each  other.  It  was  noted  earlier  that  eddy  correlation 
theory  requires  that  these  pieces  of  equipment  should  be  as  close 
together  as  possible,  ideally  sampling  at  the  same  place.  However, 
because  the  ASASP  was  only  being  sampled  at  4Hz,  this  cause  of  error  is 

unlikely  to  have  been  the  major  component  of  the  calculated  uncertainty. 

Having  said  this,  It  would  still  be  desirable  in  the  future  tc  have  a 

sample  tube  leading  from  the  centre  of  the  sonic  anemometer  to  the  ASASP 

Inlet.  This  would,  of  course  Introduce  a  time  delay  between  the  ASASP 
and  the  anemometer,  but  this  can  be  removed  when  processing  the  data. 

The  major  cause  of  error  Is  thought  to  have  been  the  low  sample 
volume  of  the  ASASP,  and  doubts  have  been  raised  as  to  the  suitability 
of  the  ASASP  to  eddy  correlation  studies.  It's  response  time  Is  slow 
compared  with  other  particle  measuring  systems  (such  as  the  ASASX,  also 
manufactured  by  PXS) ,  and  the  sonic  anemometer.  The  ASASP  Is  mere  suited 
to  measuring  dry  deposition  using  the  gradient  technique. 

To  be  able  to  draw  any  conclusions  on  the  dry  deposition  of  aerosol 
particles  would  require  further  experimentation.  In  particular,  periods 
of  more  constant  conditions  would  be  needed,  and  comparison  between 
several  days  of  data  would  be  useful. 
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Ms-  1987 


16:47:7?  to  17:4?)?5  or,  75 
Error  in  hulk  f  1  <  i  v  is  3  8  4 . 15  / 

Error  u,  hulk  deposition  velooits  is  0.705 
17:55:  6  to  18:55:40  or,  75  Mss  1987 
Error  in  hulk  flu:  is  .706.9.1 

Error  in  0,1)1-  deposition  v  e  l  o  r  1 1  s  is  0.513 
18)55:45  t.o  19:56:44  on  75  Mss  1987 
Error  l  r,  hull-  flu-,  is  3  6 .  9  3 

Error  ir,  hul)  deposition  vplorits  is  0.55  7 
19:56:49  to  ■‘0:57:33  on  75  Mss  J987 
Error  ir,  hulk  flu:  is  355.85 

Error  in  bulk  deposition  vplorits  is  0.583 


SECT  I OH  5 


KEASUREKEHTS  OK  AEROSOL  SIZE  DISTRIBUTIOHS  AT  GREAT  DUH  FELL 

a.  Intercomparison  of  Probes 

Prior  to  this  experiment  the  UMIST  Knollenberg  ASASP  was 
compared  in  detail  with  the  University  of  Galway  Knollenberg  ASAS-X 
probe.  These  experiments  were  performed  in  the  laboratory  at  UMIST  and 
involved  lntercompartsons  using  the  background  aerosol,  latex  spheres  of 
known  size  and  an  ammonium  sulphate  aerosol  distribution  generated  from 
a  nebuliser.  The  results  of  this  intercomparison  are  summarlood-tn—Iahl-e 
rmnid  will  be  used  to  compare  the  results  obtained  by  the  University  of 
Galway  with  those  obtained  by  UMIST. 

b. The  Field  Measurements 

Measurements  were  made  simultaneously  by  the  two  institutions 
during  the  period  19-20  April  1988.  The  air  trajectories  passed  close  to 
both  sites  on  20  April  when  the  transit  time  between  tf.  Ireland  and 
Great  Dun  FeLl  was  lbout  10  hours.  The  UMIST  ASASP  was  located  at  the 
site  of  our  mobile  laboratory  which  was  located  on  the  SV  face  of  Great 
Dun  Fell  approximately  150m  below  the  hill  summit  of  847m.  The  ASASP  was 
mounted  1.0m  above  the  ground  and  pointed  into  wind.  Measurements  of 
wind  speed,  wind  direction  wet  and  dry  bulb  temperature  were  made.  The 
results  obtained  are  shown  in  table  2.  The  next  step  in  analysing  this 
data  will  be  to  compare  the  results  with  those  obtained  by  the 
University  of  Galway. 
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CONCLUSIONS 

Measurements  have  been  made  oi  the  deposition  rate  oi  aerosol 
to  the  suriace  by  tne  occult,  dry  and  rain-out  processes.  It  has  oeen 
demonstrated  that  the  deposition  rates  oi  chemical  agents  when  dissolved 
in  cloud  droplets  are  large,  the  deposition  velocity  being  close  to  tnat 
for  momentum  to  the  same  suriace  when  the  cloud  is  in  contact  with  the 
ground.  These  values  are  at  least  lux  greater  than  the  dry  deposition 
velocity  of  natural  aerosol  to  the  same  surface,  fhis  result  suggests 
that  a  considerable  enhancement  of  the  deposition  rates  oi  soluble 
chemical  agents  Is  to  be  expected  in  elevated  regions  with  complex 
terrain.  The  work  has  also  shown  that  the  effeciency  of  the  rain-out 
process  is  substantially  increased  over  elevated  topography  where  the 
chemical  agents  are  dissolved  in  the  low  level  hill  cap  clouds.  This 
process  will  give  deposition  rates  oi  the  chemical  species  lux  higher 
than  the  occult  deposition  rates  although  over  the  limited  periods 
during  which  rainfall  is  actually  occurring. 


